The low-temperature collapse of the Fermi surface and phase transitions in correlated 

Fermi systems 
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A topological crossover, associated with the collapse of the Fermi surface in strongly correlated 
Fermi systems, is examined. It is demonstrated that in these systems, the temperature domain where 
standard Fermi liquid results hold dramatically narrows, because the Landau regime is replaced by 
a classical one. The impact of the collapse of the Fermi surface on pairing correlations is analyzed. 
In the domain of the Lifshitz phase diagram where the Fermi surface collapses, splitting of the BCS 
superconducting phase transition into two different ones of the same symmetry is shown to occur. 
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Introduction. For the past few years, the investiga- 
tion of topological transitions in correlated Fermi sys- 
tems that dates back to a pioneer work by I. M. Lifshitz, 
published in 1960^ has become one of hot topics in con- 
densed matter physics^— These transitions are respon- 
sible for non-Fermi-liquid (NFL) behavior that manifests 
itself in singularities of thermodynamic characteristics of 
Fermi systems. E.g. in FL theory, the spin susceptibil- 
ity x(^) remains unchanged at T — > 0. Contrariwise, at 
the transition point, \(T —> 0) diverges as T~ a , with the 
critical index a > 0, somewhat depending on the shape of 
the single-particle spectrum e(p) i 11 ' 16 i 17 In case the func- 
tion e(p) is measured from the chemical potential /z, the 
topological rearrangement of the Landau state is associ- 
ated with the change in the number of its zeroesii^ In 
conventional nonsuperfluid homogeneous Fermi systems, 
whose Fermi surfaces are singly connected, equation 



e(p) = 



(1) 



has the single real root, the Fermi momentum pp. In 
the original article, Lifshitz analyzed topological transi- 
tions, occurring in noninteracting electron gas of met- 
als at high pressures. Scenarios for such transitions, 
entailed by interactions between quasiparticles, one of 
which is addressed in this article, have emerged thirty 
years later (for a recent review, see Ref|2J). NFL 
behavior of Fermi systems near one of the critical points 
of Eq.flT]) associated with the divergence of the effective 
mass M* , (this point is called the quantum critical point 
(QCP)), has been studied extensively in recent years i 25 : 26 
Usually Eq.([T]) is analyzed at zero temperature. How- 
ever, the spectrum e(p) depends on T, being a functional 
of the quasiparticle momentum distribution n(j>, T) that 
has the standard Fermi-likc form 



n(p,T) = (l + e^ 



normalized in 3D by ordinary condition 

n3 



J n(p)dv = p =E!L ) 



(2) 



(3) 



with the volume element in 3D momentum space dv = 
p 2 dp/ir 2 . 

The key point of this article is that at T > where 
topological transitions are replaced by crossovers, there 
is a different route of the topological rearrangement, as- 
sociated with the collapse of the Fermi surface. The 
collapse occurs in case Eq.([T]) has no real roots at all. 
A classical Maxwell reconstruction of the T = ideal- 
Fermi-gas momentum distribution is the simplest exam- 
ple of such a topological rearrangement. In this case 
where e(p) = p 2 /2M — n(T), the trajectory of the sin- 
gle root of Eq.Q, denoted further px/2 due to relation 
n {Pi/2) = 1/2, stemming from Eq.@ at e(p) = 0, is eas- 
ily traced. At T — > 0, one has Pi/a(0) = pf- With the 
temperature rise, P\/2{T) moves toward the origin, at- 
taining it at temperature where /i(T) changes sign. At 
higher T, real roots of Eq.([T]) no longer exist, and the 
function e(p) becomes positive defined, the gap in the 
spectrum e(p) increasing with the further increase of T 
that entails rapid spread of the quasiparticle momentum 
distribution n(p,T). 

In systems with weak and/or moderate correlations, 
such an option is of no interest, since temperature Tm, 
at which the Maxwell-like crossover takes place, is com- 
parable with the Fermi energy = p 2 F /2M. How- 
ever, with strengthening correlations, the effective mass 
M* increases, and correspondingly, the FL spectrum 
e {p) = Pf(p — Pf)/M* becomes natter and natter. As 
a result, temperature Tm goes down that leads, in its 
turn, to the dramatic shrinkage of the temperature re- 
gion where the standard FL behaviors: C(T) oc T, for 
the specific heat, and x(T) = const, for the spin sus- 
ceptibility, hold. Importantly, at T ~ Tm -C e^ 5 the 
Landau-Migdal quasiparticle pattern remains applicable 
to evaluation of thermodynamic properties. This allows 
one to apply the quasiparticle formalism to the inves- 
tigation of relevant problems, like splitting of the BCS 
superconducting phase transition into two ones, (see be- 
low), notwithstanding discrepancies between experimen- 
tal data and standard FL predictions. 
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Simple model of the collapse of the Fermi surface at 
low T. Within the quasiparticle picture, all the ther- 
modynamic properties of correlated Fermi systems are 
evaluated in terms of the single-particle spectrum e{p). 
This spectrum can be calculated on the base of FL equa- 
tion, connecting e(p) with the FL quasiparticle momen- 
tum distribution in terms of the first harmonic f\ of 
the interaction function /, taken as a phenomenological 
input. In the 3D case, this equation has the form: 



dp 



P_ 

M 



fi{p, Pi 



dn{pi 
dpi 



■dv\ 



(4) 



The interaction function /(p,Pi) is known to coin- 
cide with a specific limit of the scattering amplitude 
r(pi,P2;q,w) where q = p — Pi — > 0,uj —> 0;q/uj —> 0. 
Generally T contains two different components. The first 
that prevails in the vicinity of second-order phase transi- 
tions changes rapidly near the Fermi surface. The second 
varies smoothly until momenta reach values much larger 
than pf- Neglecting irregular components, the scatter- 
ing amplitude T can be written in the standard form 
r oc l/(a _1 — r e q 2 /2) where a is the scattering length, 
and r e , the effective range. This form has to be supple- 
mented by an exchange term to yield 



Tt ^ 4tt / a 

T(p,Pi,q) = 77 



a/2 



M \l-ar e q 2 /2 l-ar e (p-pi+q) 2 /2 

(5) 

To facilitate the analysis we expand this expression into 
the Taylor series. Retaining only two first terms in this 
expansion, inherent in an effective mass approximation, 
where 

/( P , Pl ) = r( P , Pl) q = 0) = ^ (l - ^i(p - Pl ) 2 

'(6) 

the group velocity de(p)/dp is evaluated from Eq.((3]) in 
the closed form: 



de(p) 



dp =ii(l-^r eP ) 



(7) 



The spectrum e(p) is then calculated straightforwardly. 
In doing so an immaterial constant, associated with the 
first term in Eq.© that contains the scattering length a, 
is absorbed into the chemical potential fi to yield 



e(p) 



2M* 



M 
M 



= 1 — 2ira r e p 



(8) 



Upon accounting for a contribution er from long- 
wavelength spin fluctuations, relevant e.g. in the case 
of 3D liquid 3 He discussed below, these formulas change. 
Nevertheless, as it was shown in microscopic calculations 
of the spectrum e(p), performed in Refl27l. the effective 
mass approximation remains adequate, merely the effec- 
tive mass is modified to 



M 
~M* 



1 — 2na r e p — a 



(9) 



with a ~ 0.5. With these results, the temperature evo- 
lution of the single root P\/2{T) turns out to be identical 
to that in ideal Fermi gas. 

Temperature Tm, at which the chemical potential 
changes sign, is straightforwardly evaluated from the nor- 
malization condition <j3j> and formula (|8|) to yield 



- M 



Pf 



2M*(p)' 



(10) 



As correlations are strengthening, the ratio M* /M in- 
creases and consequently, temperature Tm goes down. 

Unfortunately in the vicinity of the QCP where the 
effective mass M* diverges, the effective mass approxi- 
mation begins to fail, because other terms of the Taylor 
expansion come into play. For illustration, let us retain 
only the next term cx (pi — P2) 4 - In this case, the group 
velocity de/dp acquires the form 



de[p) 



P . P 
-, — — Ck— — + v—^- 
dp M 



P F M 



(11) 



with v > and a, being dimensionless quantities, whose 
values are supposed to be small. As long as a keeps pos- 
itive sign, the bifurcation momentum p\, equals 0. At the 
critical point where a vanishes, one has de(p — > PF)/dp = 
Pf/M* = up F /M, so that M/M* = v < 1. Repeating 
the same manipulations that give rise to Eq. (|10p . one 
then finds Tm oc vpp/M, and therefore Tm/(-f ^ 1- 

When a changes sign, the single-particle spectrum, 
found with the aid of simple integration of Eq. (|ll[) . can 
be conveniently rewritten in the form: 



eOr) 



Pf ( 
v — — (x 



Xb) 



(12) 



where x — p 2 /p F and Xb = \a\jv. Again the effective 
mass M* turns out to be enhanced: M* /M ~ v . 

The single root p F of T = equation (TTJ persists 
as long as x;, < 1/2, otherwise, this equation acquires 
the second root p < (0) = pF(2xb — 1), and the Lan- 
dau state is rearranged. The single-particle states with 
p < < p < pf, where pf is a new Fermi momentum, 
determined by the normalization condition ([3]), remain 
filled, while states with p < p < and p > pp turn out 
to be empty. With the rise of temperature, both the 
roots of Eq. ([IJ move to meet each other at the bifurca- 
tion momentum pb = pf^/xT- Critical temperature Tm is 
straightforwardly evaluated from Eq.([3]). In the relevant 
case Xb — 1, one finds 



M 



vp 2 F /M < e° F , 



(13) 



implying again that the elevation of temperature results 
in the rapid shrinkage of the FL domain. 

Discussion. All the analytically solved models, ad- 
dressed above, have a common feature: strengthening 
interactions results in the shrinkage of the temperature 
domain where standard FL results hold. This conclu- 
sion remains valid even if flattening of the single-particle 
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spectrum e(p) takes place only for occupied states, the 
situation, inherent in many strongly correlated systems 
where the spectrum e(p) has a more complicated struc- 
ture, ( see e.g. Ref|3). The shrinkage becomes especially 
pronounced beyond the QCP ) 4 ' 22 i 23 where as mentioned 
above, the Fermi surface becomes multi-connected. For 
illustration, let Eq.((T]) have n roots at T = 0, implying 
that the Fermi surface has n sheets. Assuming the single- 
particle spectrum to be a parabolic function between two 
neighbour roots of Eq.([T]), one finds 



where e ma x = max|e(p)| inside the Fermi volume, i.e. 
for occupied states, and £qcp is the QCP Fermi energy. 
Since, as seen from Eq.([5J, the critical collapse condition 
n{p) < 1/2 is met at Tm — £max, we infer that beyond 
the QCP, critical temperature Tm rapidly falls with in- 
creasing the number n of sheets of the Fermi surface. 

Let us now briefly discuss the collapse of the Fermi 
surface in strongly correlated systems, where at T = 0, 
there exist flat bands. A different name for this disper- 
sionless portion of the single-particle spectrum e(p) is the 
fermion condensate (FC). 19 As long as the FC density is 
small, the collapse is insensitive to the presence of the 
FC. Although with further increasing the FC fraction, 
temperature Tm rapidly falls, it keeps a nonzero value, 
even if all the quasiparticles get to the FC. This conclu- 
sion holds until fx attains a critical value, at which the 
inequality n(p) < 1/2 is met for any momentum p. In 
this situation, Eq.Q has no roots, and hence, such Fermi 
systems have no Fermi surface even at T = 0. 

It is instructive to trace how the standard FL regime 
that operates in the T — > limit gives way to a different 
one at T > Tm where Eq.(U) has no roots at all. For 
illustration, let us address the spin susceptibility x(T'), 
neglecting for a while spin-spin correlations. In this case, 
x(T) is given by the standard FL formula 

X(T) cx ~ J n(p)(l - n{p))dv . (15) 

In the FL regime, overwhelming contributions to this in- 
tegral come from a domain \e[p)\ < T, adjacent to the 
Fermi surface, and after simple algebra with employing 
the FL formula dp /ale — M*/pp, we arrive at the Pauli 
result: x(T) = const cx M* . 

The collapse of the Fermi surface leads to the recon- 
struction of the FL regime. As T increases away from 
Tm, the chemical potential p becomes more and more 
negative, that triggers exponential suppression of the in- 
tegral with the n 2 (p) term on the r.h.s. of Eq . (flSl) . Even- 
tually this contribution becomes negligible, and the inte- 
gral (fTS"|) is then evaluated in the closed form with the aid 
of the normalization condition ([3]), producing the classi- 
cal Curie law: x(T) oc p/T. Accounting for the spin-spin 
correlations introduces a correction, associated with the 
Weiss temperature @w> to yield x(T) cx p/(T — @w)- 
Thus from the analysis of the spin susceptibility, we infer 



that temperature Tm, at which the Fermi surface col- 
lapses, does specify the topological crossover from the 
standard FL regime to a classical-like one. The verifica- 
tion of this conclusion on the base of experimental data 
on the specific heat C(T) encounters difficulties, requir- 
ing information on the spectrum e{p, T) at T > Tm over 
a wide momentum range that remains scarce yet. 

3D liquid 3 He is one of the systems that experiences 
the topological crossover, discussed in this article. The 
Fermi energy of this liquid is e° F ~ 5K , while the ef- 
fective mass M*(p) > 3M. At very low temperatures, 
behaviors of both C(T) and x(T) are known to be con- 
sistent with FL theory. However, already at T ~ 0.2K, 
the spin susceptibility x{T) begins to decline, and at 
T > p%/2M* ~ 1.5-K" its behavior is, indeed, obeys the 
classical Curie formula x(T) cx p/T. Meanwhile conven- 
tionally, the departure of experimental data from the FL 
result x(T) = const is attributed to damping effects. 25 
Alas, how accounting for damping may lead to the Curie 
behavior of x(^) m liquid 3 He, is hardly explicable. 

Let us now turn to strongly correlated electron sys- 
tems where the effective electron mass M* is enhanced 
as well. To clear up confusion, caused by the presence of 
the lattice field in solids, it should be pointed out that it 
results in a modification of Eq.(U]) that can be done on 
the base of FL relation, ( see e.g. Refl29f). 

,r*(p) = z r"(p) + //(p, P1 )^M^, (i6) 

connecting lo— and fc-limits of the vertex part T(p), in- 
dependently of whether the external field presents or not. 
In this equation, z stands for the renormalization factor. 
The l.h.s. of this equation is expressed in the closed form: 
zT k (p) = -MdG~ x {e = 0,p)/dp = Mde(p)/dp. As for 
the quantity zT"(p), it can be written in the Migdal 
form: 

zV(p) = e g p , (17) 

where e q is the quasiparticle effective charge 2 ^ that equals 
unity in homogeneous matter due to momentum conser- 
vation. In solids, momentum conservation breaks down, 
and therefore e q ^ 1. Additionally, in solids, the electron 
mass M is replaced by the so-called LDA electron mass 
Mlda, which cannot be evaluated with great precision 
because of the complicated structure of the lattice field. 
Therefore the effective charge e q can be absorbed into 
Mlda- With this specification, the above analysis can 
be applied to data on the magnetic susceptibility of the 
heavy-fermion alloy Yb(Rhi_ x Co x )2Si2, obtained at dif- 
ferent dopings x in Refl3Ct Experimental results evidence 
for gradual flattening out of the product Tx(T), like in 
3D liquid 3 He. However, the authors 3 ^ attribute such 
a behavior to the localization of 4/-electrons, notwith- 
standing the localization is a phase transition, rather 
than crossover, and therefore the localization scenario 
fails to explain gradual flattening of the function Tx(T). 

Impact of the collapse of the Fermi surface on pairing 
correlations. Phase transitions, associated with pairing 
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correlations, are traditionally investigated on the base 
of the BCS gap equation. In the case of s-pairing, this 
equation reads 

r tanh E ( pi *> 

A(p,T) = -J V(p,pi) 2 J A(pi,T)(foi . (18) 

Here A(p) is the gap function, E(p) = e 2 (p) + A 2 (p) is 
the Bogoliubov quasiparticle energy, and V(p,pi) is the 
zero harmonic of the pairing interaction V(p, Pi), pre- 
sented by a set of Feynman diagrams, irreducible in the 
particle-particle channel. Upon neglecting its momentum 
dependence we are led to the well-known BCS result: 

A(0) = n D e- 2 ' x , (19) 

with the dimensionless pairing constant A = 
PfM*\Vl(pf,Pf)\/k 2 and rtp,, the Debye frequency. 

Critical temperature T*, at which pairing correlations 
die out, is found from homogeneous equation 

r tanh e ( pi - T *) 

D(p) = - J V(p,pi) 2 , D(p 1 )dv 1 , (20) 

determining the location of the pole of the two-particle 
Green function at the total momentum P = 0, (Thouless 
criterion). Conventional wisdom reads that T* is criti- 
cal temperature of termination of BCS superconductiv- 
ity. Indeed, in the weak-coupling limit, Eq. (|2"0"]) is derived 
from Eg . (TTS)) , setting there A — > 0, and therefore temper- 
ature T* coincides with BCS critical temperature T c . In 
this case, T c is expressed in terms of the gap magnitude 
A(0) as 

T c = 0.57A(0) . (21) 

However, with strengthening correlations, the density of 
states, proportional to the effective mass M* , increases, 
and both the quantities A(0) and T c soar up, while Tm 
goes down. Eventually, when inequality Tm < T c — 
0.57A(0) is met, the coincidence between T* and T c 
is destroyed. Indeed, at T > Tm, the Gor'kov term 
A 2 / (e + e{p)) in the mass operator S(p, e) has no longer 
pole at the Fermi surface, or equivalently at e — > 0, due to 
the fact that e(p, T) > at any momentum p. Therefore 
at T > Tm the Gor'kov term has to be treated on equal 
footing with regular contributions to S that results in re- 
covering the Dyson-like form of the quasiparticle Green 
function G(p,e) = (e — e(p))~ 1 . In this situation, the 
Meissner effect no longer exists, and superconductivity is 
terminated. 

To clarify the structure of the pairing correlations in 
the domain T M < T < T* where, according to Eq. ([20]l . 
these correlations still persist, it is worth remembering 
the familiar situation with low-density symmetric nuclear 
matter at T = where the deuteron pole in the n-p scat- 
tering exhibits itself in full force. In this case, analyzed 
in numerous theoretical studies^ - — the pairing correc- 
tion to the chemical potential /x(0) turns out to be in 



great excess of the Fermi energy p 2 F /2M. As a result, the 
Cooper condensate transforms into a Bose condensate of 
bound states of quasiparticle pairs, the boson radius tb 
being much less than the interparticle distance ro, (the 
BCS-BEC crossover). Analogous results have been ob- 
tained in the model of bipolaronic superconductivity— 
where the strength of the electron-phonon attraction is 
supposed to be large: A > 1. With decreasing density 
p or increasing A, the chemical potential fx(0) diminishes 
and eventually becomes negative, as in the classical sit- 
uation: n(0) — > — e Q /2 where e D stands for the quasi- 
molecule binding energy. 

Above critical temperature Tb , at which the Bose con- 
densate dies out, a conglomerate of moving bound pairs 
remains that coexists with ordinary quasiparticles, the 
well-known feature of hot low-density nuclear matter. 
This coexistence manifests itself in the enhancement of 
deuteron formation in deep inelastic nuclear reactions. 
In unison with this fact, a large diamagnetic response, 
observed in high-T c superconductors, is explained, at- 
tributing it to the normal-state Landau diamagnetism of 
the bipolaron system^ 

In contrast to the case of low-density matter, in dealing 
with strongly correlated Fermi systems, the energy e D , 
determined by the Bethe-Salpeter equation 

r tan h £ (P!' T ) 

D(p) = - / V(p, Pl ) 2T D( Pl )dv u (22) 

J 2e(pi,T) - e (T) 

depends on temperature, and as stems from comparison 
of Eqs.UUD and (J22JI, e Q (T) vanishes at T = T*. 

Our first goal is to evaluate a critical point of the Lif- 
shitz phase diagram for systems with pairing correlations, 
where both the Cooper condensate and the conglomerate 
of bound pairs disappear simultaneously that occurs pro- 
vided temperatures Tm and T* coincide with each other. 
The corresponding critical value A c is found with the aid 
of the Thouless criterion (|20"]) . In doing so we employ 
the above model with the single-particle spectrum ([12")) 
where the root Pi/2(Tm) is assumed to be located quite 
far from the origin: Pi/2(Tm)/pf — 1- Remembering 
that at T = T M oc 1/M*, the chemical potential /z(T) 
vanishes, we arrive at the following equation 

, oo 

/ e p f tanh z , . 

1cc Hm^J -^rr dz > (23) 



to find that in the case T* = Tm, 

A c -0(1). (24) 

This result holds beyond the QCP as well. Indeed, 
by definition, all new pockets of the Fermi surface 
that emerge beyond the QCP disappear at T = Tm- 
Near the bifurcation momentum P\/2{Tm), the single- 
particle spectrum e(p, Tm) has the same parabolic form: 
c{p,Tm) oc w{p—p 1 /2{Tm)) 2 as in Eq. (fT2|) . with the pref- 
actor w oc 1/M*. Evidently, upon inserting these rela- 
tions into Eq. (l2"0)) the result (|2~4"|) is recovered. 
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FIG. 1: Color online. The phase diagram of a Fermi sys- 
tem with pairing correlations vs. the effective mass M* . 
The strength of pairing interaction V is supposed to be small 
enough to avoid the regime of bipolaronic superconductivity. 
The line, separated the state with pairing correlations from 
the conventional FL state, is drawn in green. The pseudogap 
phase is separated from the BCS one by the solid blue line 
where the collapse of the Fermi surface occurs. The dashed 
blue line separates the Landau states from the normal state 
without the Fermi surface. Both temperature T and the ef- 
fective mass M* are given in arbitrary units. 



Evidently, in case M* increases, temperature 7m 
drops, while T* grows. Thus in the domain of the phase 
diagram where A ~ 1, the BCS superconducting phase 
transition is split into two ones. The lower transition 
continues to be associated with termination of supercon- 
ductivity, while the upper one turns out to be related to 
the extinguishment of the bound pairs. 

The excess T* — Tm can be shown to obey a linear 
relation 



T* - T M = T M {\ - Ac) 



(25) 



determining the range of the region Tm < T < T* where 
the quasiparticles coexist with the quasimolecules, while 
the BCS correlations are completely suppressed. 

It is instructive to evaluate the binding energy e (A, T) 
at critical temperature Tm- At small difference (T* — 
Tm) << T* , this energy is calculated from a set of two 
equations, one of which is Eq. ([23|) . while the second reads 



1 cx A-i 



tanh z 



M*T M J zV 2 {z + e /T M ) 
o 



dz 



Straightforward calculations then yield 
e (X,T M ) oc 7m (A - A c ) 2 



(26) 



(27) 



One of consequences of this result is that the pairing cor- 
relations cannot regain its original BCS form at T = Tm, 
and therefore the BCS critical temperature T c must be 
lower than Tm- The analysis of the interplay between 
the two types of these correlations at low T < Tm, dat- 
ing back to Reff3ll. is not properly performed yet. I hope 
to revert to this question in the future. These results are 
summarized in the Lifshitz phase diagram of a Fermi sys- 
tem with pairing correlations, drawn in Fig. [1] where the 
effective mass M* changes, while the pairing interaction 
V remains unchanged. 

Let us now turn to pairing correlations in systems with 
flat bands, being nowadays the subject of wide spec- 
ulation because of the huge density of states of these 
systems i&iiii As we have seen, if in a system with the flat 
band, there exists the point pi/ 2 where n(p) = 1/2, such 
a system possesses the Fermi surface, and one can ex- 
pect that its superconductivity is described within BCS 
theory. In systems without the Fermi surface where in- 
equality n(p) < 1/2 is met at any point of momentum 
space, the T = ground state of the system with the 
flat band, being superfluid, is a condensate of the bound 
quasiparticle pairs. 

In conclusion, in this article, topological crossover in 
homogeneous matter, associaited with the collapse of the 
Fermi surface, is analyzed. It is shown that this collapse 
leads to the replacement of the standard FL regime by a 
classical Maxwell-like one. Temperature behavior of the 
spin susceptibility of 3D liquid 3 He is demonstrated to be 
properly explained within this scenario that implies the 
substantial extention of the temperature region where 
the quasiparticle pattern of phenomena in 3D liquid 3 He 
holds. The impact of the collapse of the Fermi surface on 
pairing correlations in a domain of the Lifshitz phase di- 
agram is examined. The topological crossover, connected 
with the collapse, is shown to be responsible for splitting 
of the BCS superconducting phase transition into two 
ones of the same symmetry. The lower one is related to 
termination of BCS superconductivity, while the upper 
one is associated with the disapperance of the conglom- 
erate of the bound pairs, a key feature of the pseudogap 
phenomenon^ 
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Ministry of Education and Science, and by Grant No. 09- 
02-01284 from the Russian Foundation for Basic Re- 
search. 
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